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ABSTRACT 



We carry out high-resolution FUSE spectroscopy of the nuclear region of 
NGC 1068. The first set of spectra was obtained with a 30" square aperture that 
collects all emission from the narrow-line region. The data reveal a strong broad 
O VI component of FWHM ~ 3500 kms" 1 and two narrow O VI AA1031/1037 
components of ~ 350 kms" 1 . The C III A977 and N III A991 emission lines in this 
spectrum can be fitted with a narrow component of FWHM ~ 1000 kms" 1 and a 
broad one of ~ 2500 kms -1 . Another set of seven spatially resolved spectra were 
made using a long slit of 1"25 x 20", at steps of ~ 1" along the axis of the emission- 
line cone. We find that (1) Major emission lines in the FUSE wavelength range 
consist of a broad and a narrow component; (2) There is a gradient in the velocity 
field for the narrow O VI component of ~ 200 kms -1 from ~ 2" southwest of 
the nucleus to ~ 4" northeast. A similar pattern is also observed with the broad 
O VI component, with a gradient of ~ 3000 kms -1 . These are consistent with 
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the HST/STIS findings and suggest a biconical structure in which the velocity 
field is mainly radial outflow; (3) A major portion of the C III and N III line flux 
is produced in the compact core. They are therefore not effective temperature 
diagnostics for the conical region; and (4) The best-fitted UV continuum suggests 
virtually no reddening, and the He II 1085/1640 ratio suggests a consistently low 
extinction factor across the cone. At ~ 2" northeast of the nucleus there is a 
region characterized by (a) a strong Lya flux, but normal C IV flux; (b) a broad 
O VI line; and (c) a significantly enhanced C III flux. 

Subject headings: galaxies: individual (NGC 1068) — galaxies: Seyfert galaxies: 
active — 



INTRODUCTION 



NGC 1068 is a prototypical Seyfert 2 galaxy. Because of its proximity (z = 0.0038) and 
bri ghtness, it has been stu died in nearly every possible detail. The polarimetric observation 
by lAntonucci et al.l (119851 ) . which reveals a Seyfert 1 spectrum in scattered light, suggests 
that the nucleus and its associated broad-line region (BLR) are obscured. This finding 
provides strong evidence for the unified theory in which yiewin g angles account for the 
differences between various active galactic nuclei (j Antonuccil 1 1 9931 ) 



The nuclear region of NGC 1068 harbors a variety of astrophysical phenomena. At the 
very center of the nucleus there is a bright compact (<0"3) region commonly referred as "the 
hot spot" . Within a few arcseconds from the n ucleus, there are several bright and compac t 
clouds that coincide with knots in the radio jets (IWilson &: Ulvestadlll983l : Evans et al.lll99ll ). 
The narrow-line region (NLR) is conical in shape toward the northeast (Fig. 1), along a 
position angle of ~ 200 deg and with an opening angle of ~ 40 deg. Beyond a 6" radius the 
surface brightness drops dramatically, and emission is dominated by two ring-like filaments 
at ~ 10" and 15" from the nucleus. 

High-spatial-resolution spectroscopy of NGC 1068 has been carried out in the optical 



(|Caganoff et al. Hl991l : lUnger et alil992l:llnglis et al 



2006) as well as in the UV (|Caganoff et al. 
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2000bl : ICecil et al.ll2002l ; iGroves et al.ll2004l ) . From approxi mately 2" southwest of the nucleus 

to 4" n ortheast, emission lines exhibit multiple components ( Cecil et alll990 : Crenshaw &: Kraemer 
2000al ): (1) major emission lines consist of narrow and broad lines; (2) broad lines are approx- 
imately 2500-4000 km s _1 wide, which may be linked to those that are found in polarized light 
and believed to be reflected light from the inner BLR, and (3) narrow lines consist of a pair 
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of red and blue components. The [O III] and [N II] line profiles suggest that the separation 
of these two components varies across the conical NLR. In addition to an overall biconical 
ionization configuration, there are compact knots whose optical spectra resemble kinemati- 



cally the associated absorption li ne systems in quasars (jCecil et al.l |2002| ; ICrenshaw. et al. 



20031 ; ICrenshaw fc Kraemerll2005l ). These line-emitting knots have blueshifted radial veloci- 
ties up to 3000 kms -1 relative to the galaxy's systemic velocity, contributing mostly to the 
emission-line flux but not the continuum. Between ~2"5 and 4"5 northeast from the nucleus, 
UV line emission is redshifted relative to the s ystemic value, a p attern that is interpreted as 
the expansion of the plasma in the radio lobe (lAxon et al.lll998l ). 

Several important emission lines in the far-UV (FUV) region between 912 and 1150 A 
are observable only with specially crafted UV instruments. During the Astro- 1 mission the 
Hopkins Ultraviolet Telescope (HUT) observed NGC 1068 with 18" and 30" apertures. The 
most striking features in the wavelengths below 1150 A are the strong C III A977 and N III 
A991 lines. The line intensity ratios of C III I(A1909)/I(A977) and N HI I(A1750)/I(A 991) are 
temperature sensitive, and the derived temperature is > 25 000 K (IKriss et al.lll992l ). higher 
than the values expected for a region producing C III and N III emission by photoionization, 
The line ratios in NGC 1068 are similar to those of the Cygnus Loop supernova remnant 
(IBlair et al.l Il99ll ). suggesting a significant contribution from shock-heating mechanisms. 
Astro-2 observations of NGC 1068 were o btained with a 12" aperture at three different 
positions. The results (IGrimes et al.l Il999l ) suggest that the emission lines observed with 
HUT likely arise in the inner nuclear region imaged with HST. However, the poor angular 
resolution of the HUT instrument (> 10") does not allow a study of the NLR in terms of 
spatial details. 

The line-emitting mechanisms in the NLR, i.e. ph otoionization from the nucle us or 
shocks produced by jets, have long been under debate. iDopita fc Sutherland! (119961 ) and 
Bicknell et al.l (119981 ) proposed that emission in the NLR may be e ntirely caused by shocks. 
Velocity splitting over 1000 kms -1 , reported by lAxon et al.l ( 19981 ). in the vicinity of some 
of the bright emission-line knots provides evidence that fast shocks exist in the NLR of 
NGC 1068. However, Seyfert galaxies host powerful nuclear sources of ionizing radiation, 
and the situation can be more complex wi t h shocks or nuclear p hotoionization prev ailing 
in different environments (1 Allen et al.l Il998l ; iMorganti et al.l Il998l ) . iMorse et al.l (119961 ) sug- 
gested that photoionizing shocks are important when a radio jet interacts with the interstellar 
medium, but n o t in t he objects where sharp, straight-edged ionization cones are observed. 
Ferguson et al.l (119951 ) argued that strong C III A977 and N III A991 emission may arise 



from fluorescence in photoionized gas if turbulent velocities exceed ~1000 kms . However, 



Grimes et al.l (119991 ) found that such a velocity would lead to extre me physical conditions 
that are inconsistent with the Astro-2 data. More recent HST data (ICrenshaw &: Kraemer 
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2000bl ; ICecil et al.ll2002l ) found that the emission-line ratios are consistent with photoioniza- 



tion instead of shock heating mechanism. 

We have carried out observations with the Far Ultraviolet Spectroscopic Explorer (FUSE) 
to study the spatial distribution of FUV emission lines. In this paper we present the results 
of both large-aperture and spatially resolved spectroscopic observations with FUSE. For the 
first time, we are able to study the position dependence of several important diagnostic lines 
in the FUV band. 



2. DATA 
2.1. FUSE data 

FUSE cove rs a wavelength ra nge between 904 and 1188 A, with a spectral resolution 



of R ~ 20000 (IMoos et al.l 120000 . It is based on a Rowland circle design and consists 
of four separate optical paths or channels. Each channel consists of a mirror coated with 
aluminum plus LiF or SiC, a focal plane assembly, which includes the spectrograph apertures, 
a diffraction grating coated with aluminum plus LiF or SiC, and a portion of FUV detector 
(two named A and B). In all, data are collected from eight traces, and the pairs of A and B 
are called channels. Our original plan was to take spectra of the nuclear region of NGC 1068 
at seven positions separated by ~ 1". Because of thermal instability of the optical structure 
in orbit, the four channels in the FUSE spectrograph are not perfectly aligned, with orbit- 
dependent drifts possibly as large as 6". Therefore, a new observing strategy was developed 
to counter the drifting effect. The high-resolution aperture (HIRS: 1"25 x 20") is used for 
spectroscopy on the finest spatial scales, to ensure that maximum resolution is maintainable 
even if the telescope imaging or pointing stability degrades below specifications. Since the 
LiFl channel is mounted on the same optical system as the FUSE fine error sensor, the 
spatial knowledge and stability is fully available for the data from this channel. 

The observations of NGC 1068 were carried out with two modes between 2001 November 
28 and December 02. The first set of data of 21951 seconds was taken with fixed pointing 
and the LWRS (low-resolution) square aperture of 30". With such a large aperture, fluxes of 
the nuclear region are all collected (Fig. 2), even with the thermally- induced mirror motions 
anticipated. The second set of observations was carried out in a scanning mode across the 
conical emission-line region, with a total exposure of 94434 seconds. A narrow HIRS slit at 
a position angle of 137 deg was used. Because of orbital constraints, we were not able to 
place the slit completely perpendicular to the conical axis. The slit actually makes an angle 
of ~ 97 deg with respect to the conical axis. 
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Over years of FUSE operation, a large amount of engineering data have been collected, 
and the drift pattern among the four channels has been well studied. After peak-up align- 
ments of channels with a bright point source at the beginning of orbital nights, the relative 
drifts exhibit a dependence on orbital time (Fig. 3). The LiF2 channel exhibits relatively 
small drifts, within one arcsecond with respect to the LiFl channel. The SiC channels exhibit 
drifts that can be represented by a constant drift rate during orbital nights. For example, 
the SiC2 channel drifts relative to the LiFl channel at a rate of approximately 4" per 30 
minutes. The starting drift at the orbital sunset is however variable, depending on the degree 
of occultation because the drift is generally in the reverse direction during orbital days. 

NGC 1068 itself is too diffuse to appear in the guide-star camera. For target acquisition, 
we used an offset guide star ~ 85" from the nucleus. We also used a FUV-bright star, Feige 
23, which is approximately four degrees away from the source, for instrumental peak-up 
alignments at the beginning of each pointing. One orbit of peak-up operation on Feige 23 
was made to establish the drift pattern due to orbital motion. During the second orbit, one 
peak-up with the MDRS aperture (4" x 20") and another with HIRS were carried out on 
Feige 23. The instrument was then pointed to the southwest part of the nucleus of NGC 
1068, and we started observations in a scanning mode. As shown in Table 1, each of the 
five pointings consists of 7-9 orbits on the target. The observation window varied orbit 
by orbit, between ~ 1200 and 2600 sec. The scanning speed in each individual orbit was 
set to complete 10" within the designated observation window. After accumulating a total 
exposure time of ~ 18000 sec, FUSE returned to Feige 23 and carried out one peak-up with 
MDRS and another with HIRS, then moved back to NGC 1068 for a new set of scans. After 
a total of five pointings, FUSE returned to Feige 23 for one last orbit of peak- up with MDRS, 
to acquire additional information on the drift pattern. 

Since the channel alignments are not perfect except at the beginning of the first orbit, 
the orbit-time dependence of channel drifts derived from the engineering data is only useful 
in terms of the relative drift speed, not for the absolute timing. 

The drifts in the SiC channels can be as large as 6", which is near the size of the ionization 
cone. The drift direction is mainly along +x on the FUSE aperture plate, i.e., along the 
conical direction of NGC 1068 at the desired FUSE aperture position angle of ~ 105 deg. 
Therefore the observation procedure was designed so that each observation starts ~ 5" from 
the nucleus in the anti-cone (southwest) direction and sweeps across the cone towards the 
northeast direction at a constant angular velocity. Fig. 4 illustrates the time sequence in one 
orbit of observation: with a constant angular speed, the FUSE HIRS slit swept across the 
nuclear region of NGC 1068. We identified the peak position of the count rate in the LiFl 
channel as the nuclear position. Based on the angular speed illustrated in Fig. 3, seven time 
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intervals were defined, and the spectral extraction in these segments yields the data A-G. The 
upper panel of Fig. 4 marks the segments for the LiFl channel, whose spatial information is 
fully available. For the other three channels, their count rates peaked at different times as a 
result of drifting. We marked their respective peak positions. The sweeping speeds in these 
channels are the sum of two terms: that of the LiFl channel plus the relative drift speed of 
each channel (Fig. 3) with respect to LiFl. Their time segments A-G are also illustrated in 
Fig. 4. The observation results in a series of segments that are dependent on the relative 
positions with respect to the nucleus. In each time sequence, one can determine the time 
when the observation slits coincide with the nucleus, as determined by the peak in the count 
rate. 

The LiFl slit positions are illustrated in Fig. 1. Seven sets of data are extracted with a 
relative shift of 1 arcsec (l'/25 at the nucleus), starting from 2" southwest from the nucleus 
(in the opposite direction of the cone), to ~ 4" northeast of it. They were named as A,B,C 
(the nucleus, see Fig. 1), D,E,F and G. The significant astigmatism of FUSE yields a spatial 
resolution of ~ 5" at best, therefore it cannot be used like a normal long-slit spectrograph 
to dissect the fine details along its slit. Not all the segment data are extracted: If a time 
segment is incomplete ("G" in Fig. 4, for example), they are not used. 

Four pairs of spectra were derived from the FUSE channels. Each pair of "A" and 
"B" data is subject to the same normalization factor. We used the O VI emission flux, 
which is common to all four channels, to normalize the spectra to the level in LiFlA. We 
first calculated the continuum level between 1050 and 1070 A for LiF2B, SiClA, SiC 2B 
and LiFlA, then subtracted it from the spectrum to enable the measurement of emission- 
line fluxes. We then calculated the O VI flux in two bins: a narrow component between 
1041.1 and 1043.6 A, and a broad component between 1027 and 1048 A. The flux ratios with 
respect to their values in the LiFlA spectrum were calculated, then their mean values were 
used as the normalization factors (Table 2). For the first set of data, taken with LWRS, we 
corrected the "worm effect" in the LiFlA, B channel, according to the procedures described in 



the FUSE Instrument and Data Handbook ( http : / / fuse. pha. j hu.edu / analysis / dh book.html [ ) . 



FUSE data reduction was carried out with pipeline version CALFUSE2.3, with specific 
timing flags to extract time-dependent (hence position-dependent) segments of spectra in 
eight traces. The FUSE spatial resolution is 1'.'5, and the spacecraft jitter, as determined from 
the FUSE engineering data, is approximately 0'/6. The total spatial resolution is therefore 
l'/6. Since the data were obtained with a l'/25 slit, there is overlap between positions. This 
can be seen by the fact that the sum of the fluxes in all positions is larger than that taken 
with one large 30" aperture. 

The normalized spectra in the four channels are binned to O.lA, then merged, weighted 
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by their signal-to-noise ratio. The merged spectra at the seven slit positions are plotted in 
Fig. 5. The wavelength bin between ~ 1077 and 1087 A falls into a gap between the LiF 
channels, and only data in the SiC2B channel exist. 



2.2. HST data 

We retrieved archival HST/STIS spectroscopic data to complement our FUSE data. 
Two sets of UV spectra are availab le: one taken at a position angle of 218 deg and four 



adjacent slit positions of 0"2 x 52" (jCecil et al.ll2002l . STIS-A hereafter), and anoth er at a 



position angle of 202 deg taken with a slit of O'.'l x 52" (ICrenshaw fc Kraemerll2000bl . STIS- 
B). Both sets (see Fig. 1) contain data taken with grating G140L and G230L, covering a 
wavelength range of ~ 1150 to 3170 A at a resolution of R ~ 1000. These STIS observations 
were made between 1988 and 2000, therefore not simultaneous with the FUSE data. For 
STIS-A, we only used three slit positions as the other does not have G140L data. The fluxes 
in the three extracted spectra are summed up. 

Since the FUSE spatial resolution (> 5" along its slit) is considerably lower than that of 
STIS, and since the position angles of these two sets of observations are nearly perpendicular, 
we must normalize the STIS data to that of FUSE. Using the continuum level around 1180 
A in these spectra did not result in satisfactory matches as the line/continuum ratios are not 
constant at these spatial scales across the nucleus. We therefore use HST WPFC2 images 
retrieved from MAST to determine the normalization factor between STIS and FUSE data. 

To extract the STIS spectra, we first smoothed the 2-D images by 30 pixels (0"67) 
with a Gaussian kernel along the spatial direction of the FUSE scans, to match the spatial 
resolution of the corresponding FUSE data. We then ran the standard STIS pipeline task 
"xld" with steps of one arcsecond (1"25 at the nucleus) that correspond to the respective 
moving FUSE slit positions, and a slit width of l'/25. To normalize the STIS fluxes, we used 
a HST/WFPC2 image of NGC 1068, taken with filter F218W, and smoothed it along the 
same direction with the same kernel size as the STIS 2-D image. We then measured the flux 
in each STIS window. To normalize the FUSE fluxes, we smoothed the WFPC2 image to the 
FUSE resolution, namely 0"67 and 5", respectively, along the direction perpendicular and 
parallel to that of the FUSE slit (Fig. 1). The flux ratios in respective extraction windows 
are used to derive the normalization factors (Table 2). To the northwest direction from the 
source, there are several bright spots whose fluxes may be picked up because of FUSE's 
low spatial resolution along its slit. We tested with several extraction windows that are 
sufficiently short not to include the fluxes from these bright spots in the smoothed WFPC2 
image. Since it is unlikely that these bright spots contribute to the redshifted UV emission 
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lines, we used the flux values measured in a short window in calculating the scaling factor. 
Any uncertainties introduced by this window selection only affect the scaling of Lya and 
C IV in position A and G. 

Since the FUSE slit positions are accurate to ~ 0'.'3, we measured the fluxes in the 
neighboring regions, to make sure that the normalization factors are accurate to at least 
20%. The height of STIS extraction windows is the same as that of FUSE (~ 1"), but 
their widths (OVl-C/2) are much smaller than FUSE (20"). Therefore considerable STIS 
normalization factors were applied. The matching STIS spectra are plotted in Fig. 6. 

We compared extracted spectra from the STIS-A and STIS-B datasets, and their line 
ratios are comparable within 25%. In the following section, we only use the data from 
STIS-A, as they cover a spatial region six times as wide (0'.'2 x 3") as that of STIS-B. 



FITTING 



Spectral analyses were carried out using IRAF task Specfit (lKrisslll994j ). For the FUSE 
spectrum taken with a large aperture, we used a pair of narrow components and one broad 
component in the O VI emission line profile. Each narrow O VI component was modeled 
as a doublet whose wavelength ratio is fixed by atomic data and whose line widths are 
identical. The O VI emission feature is heavily absorbed by interstellar absorption at ~ 
1037A, therefore we introduced several absorption components. The Ly/3 emission is not 
prominent, and it is often overwhelmed by the broad O VI emission. Only in one or two 
positions is the Ly/3 emission visible. We therefore only modeled one component with a 
fixed line width. For N III and C III emission line features, one narrow component and one 
broad component were included in fitting. We fitted a power law with no extinction to the 
continuum. The fitted results are listed in Table 3. 

For the FUSE spectra taken at positions A-G, we used one narrow component plus one 
broad component to model the O VI emission, and one component for the N III and C III 
lines. To reduce the effect of channel drifts, the FUSE spectra were fitted individually for 
each channel. The fitting results in Table 4 are primarily derived from the LiFl data, which 
are not affected by channel drifts. For wavelengths below 995 A, we used the fitting results 
from the scaled SiC2A channel, which extends to 1005 A and is therefore less susceptible than 
the SiClB channel to the detector edge effects. For He II 1085, LiF2A data were used. The 
fluxes of fitted components were normalized by the values in Table 2, and then tabulated in 
Table 4. The profiles of four major emission lines at seven different slit positions are plotted 
in Fig. 7, to show the changes in the narrow and broad components as a function of spatial 
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location. 



The UV spectra suggest a low extinction level. We first fitted the continuum in the 
wavelength windows that are free of emission and absorption lines, and the results suggest 
a power-law index of a = 0.90 ± 0.04 (fx oc X a ) and an extinction of Eb-v — 0.005 ± 0.002. 
As a comparison, the HUT spectrum between 912 and 1800 A was fitted with a power-law 
continuum of a = 0.62±0.15 and Eb-v — 0.065±0.02. If we adopt this extinction value, the 
power-law index is a = 0.85 ±0.06. The similar power indices suggest that one single power- 
law continuum can fit the far-UV spectrum of NGC 10 68 without introducing signi ficant 
reddening, as suggested from the UV/optical line ratios (ICrenshaw fc Kraemerll2000al ). 



4. Discussion 

Many previous studies have gradually unravele d the kinematical complex ities of gas in 
the nuclear region of NGC 1068. The X-ray data (IKinkhabwala et al.l 120021 ) suggest that 
emission lines are formed mainly in a photoionied plasma of a temperature aro und a few 



eV. T he spatial resolution afforded by HST leads to a picture of decelerating jet (IDas et al. 



20061 ): a biconical outflow from the nuclear region sweeps up denser, ambient clouds in the 
interstellar medium (ISM) of NGC 1068. Other possibilities include ove rlapping;, discrete 
ejection that gradually dissipates (lAxon et al.l Il998l ; ICapetti et al.lll997l ). In addition, a 
high-velocity radio jet impinges on some of the clouds. Some gas expands perpendicularly 
to the axis of the jet, and the expanding radio lobe at the end of the jet also pushes on the 
ambient ISM. Near the nucleus, kinematic components span several thousand kilometers per 
second in velocity, and the continuum hot spot visible in HST images reflects a polarized 
view of the broad lines in the active nucleus. 

Our spatially resolved FUSE observations, while not at the resolution of HST, add 
information from major emission lines shortward of the HST bandpass at high spectral 
resolution. We plot the fluxes of major UV emission lines from our observations at seven slit 
positions in Fig. 8, and list the fitted line properties in Table 4. In this section, we discuss 
the four emission lines in the FUSE spectral range, along with a comparison to emission 
lines observed in the STIS spectra. 



4.1. O VI emission 



The most prominent feature in the FUSE spectra is the O VI A A 1031, 1037 emission 
line. In the data taken with LWRS, the two narrow O VI components are of FWHM ~ 350 
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kms -1 and a separation of ~ 200 kms" 1 . In the seven spectra taken with a narrow slit, 
we resolve this blend into one narrow component with FWHM of ~ 350 kins -1 and one 
broad component. In Fig. 9 we plot the O VI profiles at the seven different slit positions. 
The high spectral resolutio n of FUSE data e nables us to compare with the results of optical 
Fabry-Perot spectroscopy ( Cecil et al.lll990l ). which reveal a narrow core of ~ 300 kms" 1 . 
According to the optical data, nearly 75% of the [N II] A6583 flux is from components of 
~ 1500 kms -1 wide. Line widths at such a scale are consistent with that derived from HST 
UV spectroscopy. The narrow-line flux in the FUSE spectra is highly concentrated (60%) 
at the compact core, suggesting that the "true" NLR probably remains unresolved, at a 
sub-arcsecond scale. While the narrow O VI emission line is not the dominant component, 
its distribution is different from its broad counterpart, while only 40% of the broad O VI 
line flux is from slit position C. 

It is surprising that the O VI emission is dominated by a component that is broader 
than those seen in polarized light: in the LWRS spectrum, more than 3/4 of the O VI flux 
is from a component with FWHM ~ 3500 kms -1 that is blue-shifted relative to the narrow 
component by ~ 500 kms" 1 . This broad component is present in all seven FUSE slit spectra 
with considerable strength. It may arise from the reflected emission from the hidden BLR, 
and /o r may be the res ult of a significant velocity dispersion in the NLR. The HST FOC 
data (lAxon et al.lll998l ) reveal that emission lines near the hot knots 2" northeast (FUSE 



slit positions D and E) are split into two veloc ity systems separated by ~ 1500 kms 1 . The 



STIS spectra discussed by iGroves et al.l (120041 ) show [O III] emission knots spanning such a 
broad velocity range in the immediate vicinity of the nucleus, but not at distances of several 
arcseconds. Since the FUSE slit collects emission from a block of regions spanning several 
arcseconds perpendicular to the conical axis, the total line emission from these regions may 
be blended into one broad component. 

The ratios of C IV/O VI may provide insight into the physical conditions of the line- 
emitting regions. The value is higher for the narrow components than their broad counter- 
parts, implying a range of the ionization parameter U ~ 0.05 — 1.0 in a typical photionization 
calculation. High values of U > 1 are c onsistent with mo dels that assume the same origin 
for the associated absorbers and BLR (IKriss et al.l 120031 ) . suggesting that the clouds that 
produce the broad emission components may be of the same origin as associated absorbers. 



4.2. Velocity field 

The narrow O VI line exhibits a systematic velocity shift from position A to G by 
approximately 220 kms" 1 . This gradient in the spatially resolved spectra explains why 



- 11 - 



there are two narrow-line compon ents in the integrated f lux fro m the large- aperture FUSE 
spectrum, where it is unresolve d. ICrenshaw &: Kraemerl (l2000al ) reported a similar velocity 



pattern in their HST spectra. iDas et al.l (120061 ) successfully modeled this as a biconical 



outflow in which radial velocity changes as a function of the distance to the central nucleus: 
the emission line knots show evidence for radial acceleration to a projected distance of 2" to 
the northeast direction, followed by deceleration up to 4". The O VI line widths also increase 
at ±2" from the nucleus, probably implying a larger dispersion in these regions. 

The broad component of O VI exhibits a qualitatively similar kinematic pattern, but at 
a larger amplitude: its centroid shifts by ~ 1500 kms -1 (Fig. 10) across the same spatial 
region. Fig. 11 shows similar trends for the narrow and broad components of C IV in the 
HST spectra which have not been explicitly noted in previous st udies: the narrow component 
follows the kinematic pattern of the optical lines modeled by IDas et al.l (120061 ): the broad 
component of C IV shows behavior similar to that of O VI in the FUSE spectra. Prior 
observations that noted this blue-shifted broad component invariably attributed it to the 
reflection of the BLR. The blue shift and width (at the position of the hot spot, and in 



i ntegr ated light) are comparable to the broad polarized B./3 line observed by lAntonucci et al. 
jl985h . 



In a scattered BLR picture, a blue shift is caused by the outflowing wind from the torus 
along our line of sight, and a line width is due to the intrinsic broad line width convolved 
with the thermal width of the hot reflecting wind. As a broad component is present in 
the FUSE spectra at all seven positions, it is possible to assume that this is reflected light 
from the hidden BLR. However, a large covering factor is needed to explain the observed 
fluxes. Assuming a covering factor of 0.1, the intrinsic flux of the broad O VI emission in 
NGC 1068 would exceed that in NGC 4151. A more reasonable explanation for the observed 
broad line widths is the large velocity dispersion between bright knots. At approximately 
2" from the core (position A and E), the FWHMs of broad components are the broadest 
at 3200 kms" 1 . These maxim a coincide with the widest spliting of velocity in bright knots 



Crenshaw fc Kraemerl (2000a|). With an intrinsic dispersion of ~800 kms^ 1 and a separation 
in velocity of ~2500 kms -1 between bright knots (from STIS results), data collected by the 
long FUSE slit would exhibit a broad component of ~3200 kms -1 , which is what observed 
at positions A and E. Extended regions of hot, photoionized gas are se en in X-ray images of 
NGC 1068 JYoung. Wilson & ShopbellboOll : iKinkhabwala et al.ll2002h that could be visible 
manifestations of this hot outflow. It is n atural to assume that these high-velocity clouds 
may be related to the associated absorbers (ICrenshaw. et al. II2003I . and references therein) in 
AGN, which are mostly blueshifted. As with the lower- velocity, lower ionization emission-line 
gas, the acceleration of this high ionization gas eventually is brought to a halt by an unknown 
deceleration mechanism, which might plausibly be interaction with the ambient ISM of NGC 



1068. The evidence for deceleration at arcsecond scale may suggest that acceleration of the 
outflow materials may take place at sub-arcsecond scales. 



4.3. C III A977 

The flux ratio of C III A1909 to 977 is extremely sensitive to temperature, and the 
value measured in N1068 suggests a high temperature that is consistent with shock-heating 



( IKriss et al.lll992l ). The FUSE data taken with a large aperture (Table 3) reveal that this 
emission line consists of a narrow and a broad component. The broad component is weak 
and hence cannot be well separated in the data segments taken with a narrow FUSE slit. 
The STIS data reveal that C III A1909 emission can be fitted with a narrow and a broad 
component of FWHM ~ 900 and 3500 kms -1 , respectively. In principal, the ratio of C III 
A1909 to 977 should be calculated only between the narrow components. 

As shown in Fig. 8, the flux of C III A977 is highly concentrated in position C and 
D. High temperatures implied by this line emission may therefore be associated with the 
compact core. In a large portion of the ionization cone, the line ratios (A1909/A977) are 
considerably higher, suggesting a lower temperature. However, the line ratios are consider- 
ably lower than that derived from the HUT data, where broad C III A 1909 component was 
included in the calculation. 



4.4. N III A991 

The line ratio of N III A 1750 to A991 is also temperature dependent, and its value 
has also been used to derive a high temperature in the ionization cone. The FUSE data 
also reveal a pair of components in the data taken with a large aperture. Unlike C III 
A977, the broad component of N III A991 , as shown in Table 3, is stronger than its narrow 
counterpart. The FUSE data taken with a narrow slit shows that nearly a half of this 
broad component is in position C. In other positions, the N III A991 emission can only be 
fitted with one component of FWHM ~ 1000 kms -1 . In combination with the trend for 
C III A977, we conclude that a bulk of flux in these temperature-sensitive emission lines is 
from the compact core of NGC 1068, and their intensity is not directly tied to the physical 
conditions in the ionization cone. 
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4.5. He II 1085 

The distribution of He II A1085 is different from other emission lines in the FUSE 
spectra: as shown in Fig. 8, the He II flux varies smoothly across the NLR region, like that 
of Lya, C IV and other lines in the STIS spectra. Since He II emission is believed to be 
produced mainly by recombination and is insensitive to the gas temperature, this pattern 
of variation may simply reflect the distribution of NLR gas across the ionization cone. The 
ratio of He II emission A1640 to A1085 is believed to be a reddening indicator. As shown 
in Fig. 12, the ratio is nearly constant (~ 5.5) between slit positions B and F. This is 
consistent with the HUT result of 5.8 ± 1.6, suggesting an insignificant level of extinction 
(Eb-v < 0.05). The low values at positions A and G may suggest that the data at these 
positions are not reliable. 

The He II I(1640)/I(1085) ratio is slightly lower than that expected from recombina- 
tion, and this may suggest that the He II A1085 may be slightly contaminated by other 
weak UV emission lines. N II A1085 is a likely candidate, and the trend of decreasing 
He II (1640)/I(1085) ratio with distance from the nucleus supports this as there is a no 



ticeable decrease in ionization state at larger distan ces from the nucleus (lAxon et al.lll998 



Kraemer. et al.ll2000at iKraemer fc Crenshaw! 12000a ). As with O VI and C IV, there appears 



to be a broad component that is considerably blueshifted. However, the relevant wavelength 
range is in a gap between the LiF channels, and only SiC2B data are available. Because of 
the low S/N ratio, the reality of such a broad component is still questionable, as no such 
counterpart is found in the corresponding HST spectra of He II A1640. 



5. Summary 

We have carried out high-spectral resolution, spatially resolved spectroscopy of the 
nuclear region of the Seyfert 2 galaxy NGC 1068. The high spectral resolution of FUSE 
data enables us to study the line profiles of O VI, N III, and C III as a function of position 
across the nuclear region. Our observations using a long slit of 1'.'25 x 20" at steps of ~ 1" 
along the axis of the emission-line cone resulted in a set of seven spatially resolved spectra 
running from ~ 2" southwest of the nucleus to ~ 4" northeast. 

The O VI profiles exhibit considerable structure: a prominent broad component with 
FWHM ~ 3500 kms -1 , and a narrow component of ~ 300 kms -1 . Both components show 
a position-dependent velocity gradient along the conical axis with a velocity shift of ~ 220 
kms -1 in the narrow component and a shift of ~ 1500 kms -1 in the broad component. Both 
patterns are consistent with radial outflow from the nuclear region. Both the continuum and 



-14- 



emission lines suggest low extinction: the UV continuum is flat to the shortest wavelengths 
in the FUSE spectrum, and the best fit is a power law with virtually no reddening. The 
He II 1640/1085 ratio also is reflective of a consistently low extinction factor across the 
emission-line cone. The majority of the C III and N III emission arises in the compact core, 
suggesting that the region with extremely high temperature is very close to the nucleus and 
remains unresolved. The line emission at ~ 2" northeast of the nucleus is strong and broad, 
characterized by (a) a strong Lya flux, but normal C IV flux; (b) a broad O VI line; and (c) 
a significantly enhanced C III flux, possibly the result of hot knots associated with shock 
heating. 

W.Z. would like to thank the FUSE team, particularly W. Blair, for their painstaking 
help in planning the observations of NGC 1068 and post-observation studies of engineering 
data. 

This work has been supported in part by NASA grant NAG-8-1527 and NAG-8-1133. 
The FUSE and HST data presented in this paper were obtained from the Multimission 
Archive at the Space Telescope Science Institute (MAST). 
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Table I. FUSE Observations of NGC 1068 



Data Set 


Start Time 


Orbit 


Integration Time 


Aperture 


Peak 


-up 


(Pointing) 


UT 




(sec) 








PI 110202000 


2001-Nov-28 09:57:00 


8 


21951 


LWRS 


MDRS 




A1390201000 


2001-Nov-29 01:23:46 


9 


19541 


HIRS 


MDRS, 


HIRS 


A1390202000 


2001-Nov-29 16:34:46 


9 


19828 


HIRS 


MDRS, 


HIRS 


A1390203000 


2001-Nov-30 08:38:00 


7 


18023 


HIRS 


MDRS, 


HIRS 


A1390204000 


2001-Nov-30 20:53:34 


9 


17962 


HIRS 


MDRS, 


HIRS 


A1390205000 


2001-Dec-01 11:52:20 


8 


19080 


HIRS 


MDRS, 


HIRS 
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Table 2. Flux Normalization Factor 



Position 


LiFIA 


LiF2A 


SiClB 


SiC2B 


STIS 


A 


1.00 


1.51 


0.86 


2.48 


2.50 


B 


1.00 


1.42 


0.97 


1.70 


2.03 


C 


1.00 


1.17 


1.03 


1.07 


1.74 


D 


1.00 


1.29 


0.96 


1.29 


2.28 


E 


1.00 


1.15 


0.78 


0.87 


3.95 


F 


1.00 


1.10 


1.28 


1.17 


4.84 


G 


1.00 


1.51 


2.75 


1.83 


6.05 
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Table 3. UV Emission Lines in Spectrum with Large LWRS Aperture 



Emission Line 


Flux a 


FWHM 


Velocity 13 




A 


(10~ 14 ergss _1 cm~ 2 ) 


(kms -1 ) 


(kms -1 ) 


C III Narrow 


977.02 


46.2 ±4.0 


720 ± 52 


28 ± 16 


C III Broad 




10.1 ±5.3 


3489 ± 1626 


1047 ± 544 


N III Narrow 


990.98 


19.8 ±5.2 


934 ± 177 


381 ±47 


N III Broad 




28.5 ±6.5 


3393 ± 315 


1371 ±431 




1025.72 


28.4 ±6.5 


3392 ±314 


1340 ± 46 


VI Narrow r 


1037.63 


41.1 ±2.5 


390 ± 31 


290 ±6 


VI Narrow b 




59.7 ±3.7 


340 ±6 


76 ±4 


VI Broad 


1034.00 


338 ±8.1 


3584 ± 184 


52 ±61 


He II Narrow 


1085.15 


45 ±2 


993 ± 44 


167 ±16 


He II Broad 




14.4 ±3.2 


4286 ± 235 


-3912 ± 354 



a With E B _ V = 0.00. 

b With respect to the systemic redshift z = 0.0038 kms 
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Table 4. UV Emission Lines at Position A-G 



Line 


A 


Flux a 

(10 14 ergs s" 1 cm~ 2 ) 


FWHM b 
(km s 


Velocity b,c Comment 
(km s 4 ) 


Position A 










C III 


977.02 


2.3 ± 1.5 


564 ± 243 


604 ± 198 


N III Narrow 


990.98 


2.0 ± 1.0 


230 ± 108 


140 ± 60 


N III Broad 




0.0 ±0.0 


25434 ± 92630 


1284 ±30138 




1025.72 


0.7 ±0.4 


878 ±0 


257 ± 192 


VI Narrow 


1037.63 


7.9 ±0.6 


655 ± 38 


376 ±6 


VI Broad 


1034.00 


7.8 ±0.8 


3131 ±488 


1192 ±207 


He II 


1085.15 


3.1 ± 1.0 


1098 ± 968 


16 ± 116 


Lya Narrow 


1215.67 


37.1 ±2.8 


1580 ± 96 


393 ± 37 


Lya Broad 




19.2 ±2.1 


5283 ± 743 


393 ±0 


C IV Narrow 


1549.50 


18.8 ±6.5 


1432 ± 95 


374 ± 96 


C IV Broad 




7.6 ± 1.4 


6178 ± 1623 


507 ± 424 


He II 


1640.46 


9.8 ±0.7 


1858 ± 134 


10 ±67 


NIII 


1750.00 


0.3 ±0.5 


1818 ±1916 


3691 ± 1758 


C III] Narrow 


1908.73 


6.8 ±1.3 


2111 ±228 


43 ±72 


C III] Broad 




9.4 ± 1.3 


6883 ±1106 


-535 ± 348 


Position B 










C III 


977.02 


14.1 ±3.5 


596 ± 290 


-31 ±55 


N III Narrow 


990.98 


8.8 ±3.5 


1048 ±351 


519 ± 190 


N III Broad 




2.4 ±2.7 


2005 ± 2083 


2332 ± 1454 


Ly/3 _ ^ 


1025.72 


5.2 ± 1.0 


878 ±0 


-2 ± 105 


VI Narrow 


1037.63 


23.9 ±5.1 


385 ± 34 


333 ± 12 


VI Broad 


1034.00 


123.5 ±2.5 


2535 ± 35 


298 ± 49 


He II 


1085.15 


10.2 ±1.5 


1108 ±177 


-18 ±74 


Lya Narrow 


1215.67 


158.2 ±0.4 


1528 ± 18 


366 ± 10 


Lya Broad 




119.9 ±7.8 


7676 ± 787 


366 ±0 


C IV Narrow 


1549.50 


95.8 ±2.9 


1348 ± 32 


284 ± 12 


C IV Broad 




60.3 ±3.1 


4424 ±172 


-92 ± 66 


He II 


1640.46 


47.0 ±1.6 


1872 ± 66 


-234 ±31 


NIII 


1750.00 


4.1 ± 1.0 


1848 ± 569 


299 ± 169 


C III] Narrow 


1908.73 


29.6 ±2.4 


2229 ± 95 


-126 ±27 


C III] Broad 




36.7 ±2.2 


6296 ± 408 


-678 ±110 


Position C 










C III 


977.02 


38.6 ±4.3 


911 ±99 


-7 ±49 


N III Narrow 


990.98 


11.0 ± 10.3 


989 ± 558 


-8 ±202 


N III Broad 




29.0 ±14.8 


3207 ±767 


870 ± 770 
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Table 4 — Continued 



Line 




Flux a 


FWHM b 


Velocity b ' c Comment 




A 


(10~ 14 ergs s _1 cm -2 ) 


(km s _1 ) 


(km s _1 ) 




1025.72 


12.6 ±2.9 


878 ± 132 


-84 ± 61 


VI Narrow 


1037.63 


46.2 ±6.3 


333 ± 22 


290 ±9 


VI Broad 


1034.00 


325.3 ±43.4 


2523 ± 118 


162 ±55 


He II 


1085.15 


17.0 ±1.7 


1149 ±116 


-70 ± 55 


Lya Narrow 


1215.67 


278.5 ±4.7 


1352 ±21 


395 ±7 


Lya Broad 




253.3 ±13.5 


9104 ± 641 


395 ±0 


C IV Narrow 


1549.50 


167.6 ±9.9 


1472 ± 59 


172 ± 19 


C IV Broad 




131.0 ±9.8 


4193 ± 244 


-390 ± 73 


He II 


1640.46 


92.0 ±3.8 


1830 ± 82 


-406 ± 35 


NIII 


1750.00 


10.1 ± 1.6 


2605 ± 456 


118 ±160 


C III] Narrow 


1908.73 


38.4 ±4.4 


1902 ± 107 


-270 ± 28 


C III] Broad 




80.3 ±4.0 


5053 ±219 


-672 ± 56 


Position D 










C III 


977.02 


31.2 ±3.2 


900 ±0 


183 ±52 


N III Narrow 


990.98 


11.1 ±3.3 


888 ± 243 


372 ± 93 


N III Broad 




44.1 ± 71.1 


5775 ± 3455 


-3804 ± 4077 


Ly/3 


1025.72 


6.2 ± 1.8 


878 ±0 


-46 ±0 


VI Narrow 


1037.63 


14.9 ±0.9 


380 ±31 


287 ±9 


VI Broad 


1034.00 


181.8 ± 17.1 


2506 ± 126 


-283 ± 32 


He II 


1085.15 


15.6 ±1.2 


940 ± 85 


21 ±36 


Lya Narrow 


1215.67 


291.4 ± 14.3 


1248 ± 42 


388 ±7 


Lya Broad 




145.9 ± 10.6 


5718 ± 1708 


388 ±0 


C IV Narrow 


1549.50 


143.7 ± 18.7 


1777 ± 125 


6 ±21 


C IV Broad 




105.9 ± 18.4 


3388 ± 130 


-793 ± 17 


He II 


1640.46 


88.6 ±3.4 


1928 ± 73 


-504 ± 35 


NIII 


1750.00 


9.1 ±2.3 


2415 ± 404 


15 ± 188 


C III] Narrow 


1908.73 


45.3 ±2.9 


2181 ±493 


-367 ± 114 


C III] Broad 




58.8 ±0.0 


5201 ± 


-877 ±39538 


Position E 










C III 


977.02 


17.7 ±4.2 


1149 ±412 


27 ± 119 


N III Narrow 


990.98 


6.7 ±2.6 


577 ± 177 


375 ± 84 


N III Broad 




12.9 ± 15.0 


2584 ± 2236 


-764 ±1615 


Ly/3 


1025.72 


3.6 ± 1.2 


878 ±0 


68 ±0 


VI Narrow 


1037.63 


13.0 ± 1.2 


639 ± 42 


252 ± 17 


VI Broad 


1034.00 


52.9 ±2.9 


3145 ± 168 


-475 ± 78 


He II 


1085.15 


12.1 ±0.8 


818 ±58 


208 ± 25 


Lya Narrow 


1215.67 


359.0 ±9.4 


1265 ± 24 


474 ±7 
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Table 4 — Continued 



Line 




Flux a 


FWHM b 


Velocity b ' c Comment 




A 


(10~ 14 ergs s _1 cm -2 ) 


(km s _1 ) 


(km s _1 ) 


Lja Broad 




72.4 ±6.4 


4953 ± 1166 


474 ±0 


C IV Narrow 


1549.50 


107.2 ± 10.9 


1812 ±141 


79 ±71 


C IV Broad 




63.5 ±11.0 


3241 ± 134 


-981 ±99 


He II 


1640.46 


64.0 ±2.9 


4490 ± 157 


1112 ±5 


NIII 


1750.00 


9.1 ± 1.6 


2986 ± 549 


362 ± 237 


C III] Narrow 


1908.73 


30.3 ±4.2 


2185 ± 135 


82 ±55 


C III] Broad 




39.3 ±4.3 


4489 ± 251 


-1026 ±162 


Position F 










C III 


977.02 


15.8 ±4.3 


900 ± 270 


424 ± 119 


N III Narrow 


990.98 


7.1 ±5.8 


1056 ± 625 


164 ± 289 


N III Broad 




8.0 ±9.4 


3540 ± 1540 


333 ± 1681 




1025.72 


9.0 ±0.7 


878 ±0 


257 ±47 


VI Narrow 


1037.63 


11.4 ± 1.9 


219 ±30 


152 ± 17 


VI Broad 


1034.00 


96.1 ± 15.6 


1785 ± 16 


456 ± 20 


He II 


1085.15 


8.3 ±0.7 


665 ± 66 


266 ± 25 


Lja Narrow 


1215.67 


299.6 ±9.6 


1223 ± 33 


520 ± 10 


Lja Broad 




46.8 ± 7.4 


4582 ± 1014 


520 ±0 


C IV Narrow 


1549.50 


73.3 ±5.9 


1181 ±72 


313 ±29 


C IV Broad 




34.7 ±5.3 


3365 ± 323 


-461 ± 178 


He II 


1640.46 


42.6 ±2.8 


1132 ± 74 


305 ± 33 


NIII 


1750.00 


4.7 ± 1.6 


1289 ± 702 


1040 ±211 


C III] Narrow 


1908.73 


29.1 ± 1.2 


1500 ± 77 


400 ± 41 


C III] Broad 




14.3 ±0.9 


3848 ± 392 


-1679 ±139 


Position G 










C III 


977.02 


12.5 ±5.4 


344 ± 133 


91 ±80 


N III Narrow 


990.98 


8.7 ±4.2 


833 ± 398 


417 ± 148 


N III Broad 




10.8 ±58.4 


3237 ± 7653 


-2085 ± 9897 




1025.72 


4.3 ±0.9 


878 ±0 


-61 ±93 


VI Narrow 


1037.63 


13.3 ±2.6 


372 ± 16 


189 ± 14 


VI Broad 


1034.00 


37.6 ±4.9 


2852 ± 422 


416 ± 95 


He II 


1085.15 


7.5 ±1.4 


617 ± 105 


150 ± 58 


Lja Narrow 


1215.67 


152.2 ±5.8 


1208 ± 35 


523 ± 12 


Lja Broad 




32.2 ±3.9 


4627 ± 848 


523 ±0 


C IV Narrow 


1549.50 


50.8 ±3.0 


1227 ±62 


251 ±27 


C IV Broad 




19.5 ±2.3 


4877 ±581 


-283 ±216 


He II 


1640.46 


28.0 ± 1.7 


1047 ± 58 


265 ± 29 


NIII 


1750.00 


1.5 ±0.8 


708 ± 334 


1453 ± 195 
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Table 4 — Continued 



Line 


A 


(io- 


Flux a 

-14 er ggg-l Cm - 2 ) 


FWHM b 

(kms" 1 ) 


Velocity b,c Comment 
(kms^ 1 ) 


C III] Narrow 
C III] Broad 


1908.73 




12.5 ±3.1 
12.4 ±3.5 


1797 ±252 
4491 ± 658 


333 ± 83 
-870 ±510 


a With E B _ V 


= 0.00. 










b Values with 


zero errors 


are 


pre- fixed. 







c With respect to the systemic redshift z = 0.0038 kms 
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Figure Captions 



Fig. 1.— HST image of the central region of NGC 1068 obtained with WFPC2/F218W, with 
a field of view of ~ 25". The north is up, east is to the left, and the nucleus is marked with 
a cross. Seven FUSE slit positions are marked with letters A-G, while C is at the nucleus. 
Note that the actual slit width is slightly wider than that in the illustration. The narrow 
slits along the near-perpendicular direction are that of HST/STIS: The three paralle l slits 



along the northeast and southwest direction are (X'2 wi de fSTIS-A. ICecil et a 



the single long slit is another that is 0"1 wide (STIS-B, I Crenshaw & Kraemer 



2000bf ) 



2002), and 



Fig. 2. — FUSE spectrum of NGC 1068 at full aperture (30" square). The data are binned 
to 0.1 A. Propagation errors are plotted in the lower panel. The effect of major geocoronal 
emission lines are marked with the Earth symbols, and data at these wavelength bins are 
removed. 

Fig. 3. — Typical relative drifts of the FUSE SiC2 channel vs. orbital time (after orbital 
sunset), with respect to the LiFl channel. Data are collected from various orbits with 
occultation. 

Fig. 4. — Count rate vs. orbital time in a sample data set A1390201001. Different channels 
register their peaks at different times, which are presumed to coincide with the nucleus 
position. Incomplete segments are not used (Position A in SiC channels and position G in 
LiF channels). 

Fig. 5. — Merged FUSE spectra at seven FUSE slit positions, binned to 0.5 A. The nucleus 
is at position C. 

Fig. 6. — STIS spectra extracted at different windows that correspond to the respective 
FUSE slits. Fluxes are normalized to the corresponding FUSE data. 

Fig. 7. — Line profiles at different FUSE slit positions. The lines are redshifted O VI A1037 
(centered at kins" 1 ), C IV A1549, N III A991, and C III A977, respectively. 

Fig. 8. — Intensities of major emission lines at different FUSE positions, from STIS (left 
panels) and FUSE (right panels). When a line is fitted with multiple components with 
reasonably accuracy, the flux presented is their sum. For C III A977, only one narrow 
component is fitted. 



Fig. 9. — Emission and absorption lines near the O VI wavelength, from position A to G 
after subtraction of the fitted continuum. The spectrum marked "All" in the lower bottom 



-25 - 



panel is from the data with a large LWRS aperture. The dashed lines mark the wavelengths 
of redshifted Ly/3, O VI AA1031.95/1037.63, respectively, from left to right. 

Fig. 10. — Properties of fitted O VI emission line at different FUSE slit positions. Line fluxes 
are in units of 10~ 14 ergss" 1 cm -2 , FWHM and velocity are in units of kms -1 . 

Fig. 11. — Properties of fitted C IV emission line (from HST/STIS) at different FUSE 
positions. Line fluxes are in units of 10 -14 ergs s -1 cm -2 , FWHM and velocity are in units 
of kms -1 . 

Fig. 12. — Line ratios at different slit positions. The fitting results of C IV and O VI are for 
their narrow-line components. 
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